Magnetic behavior of nanoparticles in patterned thin films 
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The magnetic behavior of truncated conical nanoparticles in patterned thin films is investigated 
as a function of their size and shape. Using a scaling technique, phase diagrams giving the relative 
stability of characteristic internal magnetic structures of the particles are obtained. The role of 
the uniaxial anisotropy in determining the magnetic properties of such systems is discussed, and a 
simple method for stablishing its strength is proposed. 



Nanolithography techniques combined with material 
deposition and pattern transfer processes have made pos- 
sible the production of regular arrays of magnetic parti- 
cles with individual linear dimensions in the submicron 
range. A relatively high degree of precision has been 
achieved in the growth process of these structures, allow- 
ing the control of both the particle size and the spacing 
between them. The samples thus produced exhibit rather 
interesting properties , an d show great potential for tech- 
nological application*^. 

Electrodeposition, and evaporation followed by lift-off 
are the two methods that have been most widely used 
in the preparation of this sort of arrays or patterned 
thin films'. The first is suitable for producing cylindri- 
cally shaped particles in which the ratio R between their 
height H and diameter D of the bottom basis (aspect 
ratio) is bigger than 1. On the other hand, evaporation 
and lift-off process is more suitable for the production 
of periodic arrays covering large areas, which may reach 
several square Taking under consideration its low 

cost of production, this method is seen as quite conve- 
nient for commercial applications. Particles pr odu ced by 
evaporation have the shape of truncated cones^^, which 
can be characterized by three parameters, namely H, D, 
and the ratio £ = (D — d)/2H, where d is the diameter 
of the top basis. 

It is generally accepted that particles of ferromagnetic 
materials with linear dimensions of the order of 10 2 nm 
or less are too small to accommodate a domain wall. 
As a consequence, one could envisage the production of 
arrays of single domain ferromagnetic particles, which 
could be used in the design of spintronic devices or as 
magnetic media for high density recording^! In recent 
years, several groups have concentrated their attention 
on the magnetic behavior of nanosized particles^Hin] It 
has been found that the internal magnetic structure of 
sufficiently small particles is indeed rather close to that 
of a fcrromagnet monodomain. Larger particles, how- 
ever, may also exhibit vortex-like configurations, depend- 
ing on their size and shape. Clearly, the determination of 
the conditions under which these distinct magnetic con- 



figurations or phases are found is of great relevance for 
practical applications. However, from the experimental 
point of view, this is not a simple task since transitions 
between magnetic configurations as t he shape or size of 
the particle changes are not sharpen Hence, theoretical 
investigations of the magnetic behavior of such nanopar- 
ticles become highly desirable. 

In this paper, we focus our attention on patterned thin 
films produced by evaporation followed by lift-off. The 
systems we have in mind are Ni films, wh ich have been ex- 
tensively studied in recent years^*- 11 * 12 *. We aim at deter- 
mining, for various shapes and sizes of the particles, the 
configuration of lowest energy out of three characteristic 
ones. These configurations are: (I) ferromagnetic with 
the magnetization parallel to the cone basis (in-plane); 
(II) ferromagnetic with the magnetization perpendicular 
to the cone basis (out-of-plane) ; (III) vortex-state with 
the magnetic moments laying parallel to the cone basis. 
We shall deal with the situation in which the distance 
L between the particles in arrays is large enough (i.e. 
L > D) for the interaction between them to be safely 
neglectecP. In such cases, the magnetic structure of each 
particle is determined solely by the internal interactions 
its magnetic moments are subjected to, namely the short- 
ranged exchange interaction, the classical dipolar inter- 
action (which is responsible for the shape anisotropy), 
and the magneto-crystalline anisotropy. It is worth men- 
tioning that transmission electron microscopy has clearly 
established that both Ni and Co particles in samples pre- 
pared by evaporation are polycrystalline, with elongated 
columnar structures^. Micromagnetic calculations of the 
remanence in those particles indicate that the inclusion 
of an out-of-plane crystalline anisotropy, in addition to 
the magnetocrystalline anisotropy, is necessary in order 
to reproduce the experimental data~ The importance of 
the uniaxial anisotropy in stabilizing configuration II is 
discussed towards the end of this article. 

Thus, we consider truncated conical Ni nanoparticles 
consisting of N magnetic moments occupying sites of an 
underlaying fee structure. We assume that the cone axis 
(growth direction) is in the (111) directiorP. For fixed 
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C, the phase diagram giving the relative stability of the 
three configurations described above is obtained by com- 
paring the corresponding total energies at each point in 
the DxH plane. For each configuration of the mag- 

netic moments, it is given by 

£ tot ({7^}, = \j^[Eij-Jijlti''tfj] + UK, (1) 

where Eij is the dipolar interaction between magnetic 
moments ~~jti and ~ftj, Jy is the exchange coupling, and 
Uk is the crystalline uniaxial anisotropy term given by 
Uk = —KJ2 i=1 cos 2 9i. Here, 0% is the angle between 
~jti and the cone axis. For the systems under consid- 
eration, we take the lattice parameter clq = 3.52 A, 
Hi = H = 0.62/ig. and Jij = J = 0.9 mRy for first 
nearest-neighbors^ and zero otherwise. Estimates of the 
anisotropy constant K fall in the range 1.8xl0 5 -2.6xl0 5 
erg/cm 3 [3,11], so we take for the moment the value 
2 x 10 5 [3]. 

The calculation of the above expression is a well de- 
fined problem. However, with the present day stan- 
dard computational facilities, it becomes rather time- 
consuming as soon as the number N of magnetic mo- 
ments assumes values comparable to those in real sys- 
tems, i.e. 10 s or even 10 9 . As a consequence, the com- 
putational effort necessary for obtaining a phase diagram 
in which particles in such size range are included turns 
out to be prohibitively large. Nevertheless, it can be 
significantly reduced with the use of a scaling technique 
recently proposed and used to study the magnetic prop- 
erties of cylindrically shaped particles^. Such technique 
involves two steps: (a) calculations are performed for sys- 
tems with the exchange interaction J scaled down by a 
factor x < 1; (b) the H and D axes in the resulting 
phase diagrams are scaled up a factor (1/x)'', where 77 is 
a positive constant (see below). 

In which follows, we first demonstrate that the scaling 
technique can be also used to obtain the phase diagram 
of truncated conical particles, even in the presence of a 
strong uniaxial anisotropy, and show how 77 is obtained. 
Having done this, we use such technique to obtain phase 
diagrams for different values of £. Results are compared 
to existing experimental data. Finally, we comment on 
the effects of the uniaxial anisotropy on the magnetic 
behavior of those particles. 

Fig.(l) shows results for cones with £ = 0.3, and scal- 
ing factors x — 0.010 (squares), 0.015 (circles), and 0.020 
(triangles). Dashed line corresponds to the maximum 
height Umax = D/2C, one particle can reach for each di- 
ameter D. As in the case of cylinders^, the curves in the 
three phase diagrams have similar shapes, with a triple 
point that moves along a common line as x is varied. 
From the plotting of the height H t at the triple point as 
a function of x, we find that these two quantities are re- 
lated by the equation H t = Ax v , with A « 74.4 nm and 
77 « 0.55. The uncertainty in the values of both A and 
77 is smaller than 2%. Within the numerical accuracy of 



the present calculations, the value of 77 that we obtained 
in this case coincides with that for cylindrically shaped 
particles^. It means that 77 is either independent of £ or 
has a rather weak dependence on such parameter. This 
point has been confirmed by the calculations we have 
performed for other values of C- Having determined the 
value of 77, we then scale up the three diagrams by divid- 
ing the D and H axes by the corresponding x to power 
77. Results are shown in Fig. (2), where all points fall on 
a single set of lines, representing the phase diagram for 
the full strength of the exchange interaction (x = 1). 
We point out that the slope of the line separating phases 
I and II gives the critical value (R c ) of the aspect ratio 
R = H/D, which depends on the strength of the uniaxial 
anisotropy. We shall return to this point below. 
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FIG. 1. Phase diagrams giving the relative stability of the in- 
plane ferromagnetic (I), out-of-plane ferromagnetic (II), and 
vortex-like (III) configurations of truncated conical Ni parti- 
cles with £ = 0.3, for scaling factor x = 0.01 (squares), 0.015 
(circles) and 0.02 (triangles) . Dashed line gives the maximum 
height H m ax of the particles as a function of the diameter D 
of their basis, for the specified value of £■ The inset illustrates 
the geometry of the particle under consideration. 

We are now in position to examine the phase diagrams 
for different values of C and full strength of the exchange 
coupling. Fig. (3) shows results for £ = (squares), 0.2 
(triangles), and 0.4 (circles). Dashed and dotted lines 
correspond to H max (D) for £ equals to 0.2 and 0.4, re- 
spectively. We notice that as £ increases from zero, that 
is to say, as the shape of the particles deviates from 
that of a cylinder, changes in the lines separating the 
three phases become less pronounced. In fact, the lines 
for £= 0.3 would appear in between those for 0.2 and 
0.4, and have been eliminated from the figure for the 
sake of clarity. These results are confirmed by the ex- 
perimental data of Ross et aiP for Ni samples. These 
authors have investigated the magnetic behavior of pat- 
terned thin films consisting of truncated conical parti- 
cles with 75 A < D < 122 A, 0.43 < R < 1.21, and 
0.14 < C < 0.35. They concluded that, independently of 
the value of £, transitions from configuration I to II oc- 
curs at a critical aspect ratio R c — 0.65, which is in good 
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FIG. 2. Phase diagram obtained from those in Fig.(l), with 
the corresponding axes scaled by a factor l/x v , where rj = 
0.55. Dashed line corresponds to H ma x(D). 



find that for the whole range of values of £ considered 
here, the relation between R c and K is not much sensi- 
tive to that parameter. Indeed, for £ > 0.2 and R c < 0.7, 
the curves can be regarded as independent of (. Apart 
from confirming the intuitive idea that a strong uniaxial 
anisotropy is necessary for stabilizing the out-of-plane 
configuration in particles with large £, Fig. (4) provides 
us with an extremely simple and practical method for 
determining the value of K in truncated conical parti- 
cles. The reason is the fact that the computational effort 
necessary for obtaining such curve is minimum. In fact, 
since the border line between phases I and II, for fixed 
values of K and £, is a straight line, its slope (that is to 
say, R c ) can be obtained from calculations for small par- 
ticles. Thus, as a general procedure, the strength K of 
the uniaxial anisotropy necessary for obtaining the whole 
phase diagram can be immediately determined from a 
similar graph and the measured value of critical aspect 
ratio R c , without much computational effort. 



agreement with the phase diagrams in Fig. (3). Indeed, 
according to our results, for £ > 0, transitions between 
these two phases occur for R c « 0.61. 
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FIG. 4. The critical aspect ratio R c plotted as a function of 
the uniaxial anisotropy constant K, for cylindrically-shaped 
(open squares) and cone-shaped Ni particles with £ = 0.2 
(triangles) and 0.4 (circles). 



FIG. 3. Phase diagrams for cylindrically-shaped (squares) 
and cone-shaped Ni particles with ( = 0.2 (triangles) and 0.4 
(circles). Dashed and dotted lines correspond to H max (D) for 
£ equals to 0.2 and 0.4, respectively. 

Although the boundaries between phases in the di- 
agrams in Fig. (3) do not change much as £ is varied, 
the extent of region II (corresponding to the out-of-plane 
configuration) is reduced as £ increases. Such region is 
bounded from above by the line H max = D/2( , which 
means that, for a fixed value of the uniaxial anisotropy 
constant K, a critical value £ c = 1/2R C exists, above 
which phase II cannot be observed. However, we expect 
the uniaxial anisotropy to favor the out-of-plane config- 
uration. Thus, R c should decrease as K increases. This 
is confirmed by the curves in Fig. (4), where R c is plot- 
ted as a function of K for different values of £, namely 
zero (squares), 0.2 (triangles), and 0.4 (full circles). We 



In conclusion, we have investigated the magnetic be- 
havior of truncated conical particles, which exhibit a 
strong uniaxial anisotropy along the cone axis. We have 
shown that even in such cases, the phase diagrams giving 
the relative stability of three characteristic arrangements 
of the magnetic moments within each particle can be ob- 
tained using a scaling method recently proposed. We 
have found that the stability of the ferromagnetic out-of- 
plane configuration strongly depends on the strength K 
of the uniaxial anisotropy, particularly when the shape of 
the particles deviates significantly from that of a cylinder 
(i.e. for large £). On the basis of our results, we have 
proposed a simple method for determining the value of 
K from the measured critical aspect ratio R c . We believe 
the present work represents a relevant contribution to the 
understanding of the magnetic properties of nanoparti- 
cles, and expect it to help the development of new mag- 
netic devices. 
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